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Abstract 

In this paper we show that the dissipative version of the laser model proposed by Alii 
and Sewell can be obtained by considering the stochastic limit of the (open system) 
hamiltonian introduced by Hepp and Lieb in their seminal work. We also prove that 
the Dicke-Haken-Lax hamiltonian produces, after the stochastic limit is considered, the 
generator of a semigroup with equations of motion very similar to those of Alli-Sewell, 
and coinciding with these under suitable conditions. 



PACS Numbers: 02.90.+p, 03.65.Db 



I Introduction 



In two recent papers, [1, 2], a dissipative laser model has been introduced and analyzed in 
some details. In particular in [1] (AS in the following) the rigorous definition of the unbounded 
generator of the model, which consists of a sum of a free radiation and a free matter generator 
plus a matter-radiation term, is given and the existence of the thermodynamical limit of the 
dynamics of some macroscopic observables is deduced. Moreover, the analysis of this dynamics 
shows that two phase transitions occur in the model, depending on the value of a certain 
pumping strenght. In [2] the analysis has been continued paying particular attention to the 
existence of the dynamics of the microscopic observables, which are only the ones of the matter 
since, in the thermodynamical limit, we proved that the field of the radiation becames classical. 
Also, the existence of a transient has been proved and an entropy principle has been deduced. 

On the other hand, in a series of papers [3, 4] culminating with the fundamental work by 
Hepp and Lieb [5] (HL in the following) many conservative models of matter interacting with 
radiation were proposed. In particular, in [5] the authors have introduced a model of an open 
system of matter and of a single mode of radiation interacting among them and with their 
(bosonic) reservoirs, but, to simplify the treatment, they have considered a simplified version 
in which the matter bosonic reservoir is replaced by a fermionic one. In this way they avoid 
dealing with unbounded operators. This is what they call the Dicke-Haken-Lax model (DHL 
model in the following). 

In [1, 2] the relation between the AS model and a many mode version of the HL model 
is claimed: of course, since no reservoir appear in the semigroup formulation as given by [1], 
this claim is reasonable but it is not clear the explicit way in which HL should be related to 
AS. In this paper we will prove that the relation between the two models is provided by (a 
slightly modified version of) the stochastic limit (SL), [6] and reference therein. In particular, 
if we start with the physical AS system (radiation and matter) and we introduce in a natural 
way two reservoirs (one is not enough!) for the matter and another reservoir for the radiation, 
then the SL of the hamiltonian for the new system constructed in this way returns back the 
original AS generator, under very reasonable hypotheses. Moreover, the model which we have 
constructed ad hoc to get this generator surprisingly coincides with the HL laser model, [5]. 
This is the content of Section 3, which follows a section where we introduce all the models we 
will deal with, to keep the paper self-contained. 

In Section 4 we will consider the SL of the fermionic version of the HL model, known as 
the DHL model, [5, 7]. We will find that, even if the form of the generator apparently differs 
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from the one by AS, under certain conditions on the coefficients which define the model, the 
equations of motion for the observables of the matter-radiation system coincide with the ones 
given in AS. 

In Section 5 we give our conclusions while the Appendix is devoted to summarize few results 
on SL which are used everywhere in this paper. 

Before concluding this section we wish to remark that we call the procedure proposed here 
stochastic limit even if a minor difference exists between the original approach, [6], and the one 
we will use here, namely the appearance of different powers of the over-all coupling constant A 
which appear in our hamiltonian operators. The final remark concerns our notation which we 
try to keep as simple as possible by neglecting the symbol of tensor product (almost) everywhere 
in the paper. 

II The Physicals Models 

In this section we will discuss the main characteristics of the three physical models which will 
be considered in this paper. In particular, we will only give the definition of the hamiltonians 
for the HL and the DHL models and the expression of the generator for the AS model, without 
even mentioning mathematical details like, for instance, those related to the domain problem 
intrinsic with all these models due to the presence of bosonic operators. We refer to the original 
papers for these and further details which are not relevant in this work. 
We begin with the AS model. 

This model is a dissipative quantum system, H^ N \ consisting of a chain of 2iV + 1 identical 
two-level atoms interacting with an n— mode radiation field, n fixed and finite. We build the 
model from its constituent parts starting with the single atom. 

This is assumed to be a two-state atom or spin, S at . Its algebra of observables, A a t, is that 
of the two-by-two matrices, spanned by the Pauli matrices (cr x , o- yi a z ) and the identity, /. They 
satisfy the relations 

a l = a l = °l = J ; a x a y = io- z , etc. (2.1) 
We define the spin raising and lowering operators 

<7± = ^(a x ±ia y ). (2.2) 

We assume that the atom is coupled to a pump and a sink, and that its dynamics is given by a 
one-parameter semigroup {T at (t)\tE~R + } of completely positive, identity preserving contractions 
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of Act, whose generator, L at , is of the following form. 

L at cr± = -(7i=F?e)o"±; L at a z = -72 (tr* - r}I), (2.3) 

where e(> 0) is the energy difference between the ground and excited states of the atom, and 
the 7's and r\ are constants whose values are determined by the atomic coupling to the energy 
source and sink, and are subject to the restrictions that 

< 7 2 <2 7l ; -1<?7<1. (2.4) 



The matter consists of 2N+1 non- interacting copies of S at , located at the sites r = —N, . ., N 
of the one-dimensional lattice Z. Thus, at each site, r, there is a copy, E r , of S at , whose algebra 
of observables, A r , and dynamical semigroup, T r , are isomorphic with A at and T at , respectively. 
We denote by <7 r ,u the copy of a u at r, for u — x,y, z, ±. 

We define the algebra of observables, A^ N \ and the dynamical semigroup, T^l, of the matter 
to be ^)^L_ N A r and (8)^L_ N T r , respectively. Thus, A^ is the algebra of linear transformations 
of C 4JV+2 . We identify elements A r of A r with those of given by their tensor products 
with the identity operators attached to the remaining sites. Under this identification, the 
commutant, A' r , of A r is the tensor product ® s ^rA s . The same identification will be implicitly 
assumed for the other models. 

It follows from these specifications that the generator, L^ t , of is given by the formula 

where I N = {-N, -1, 0, 1, .., A^}. Here 

L r cr r ,± = -(liTic)vr,±; L r a r)Z = -7 2 (o" r , z - r)I)\ 
and L r (A r A' r ) = (L r A r )A' r \/A r eA r , A' r eA' r (2.6) 

We assume, furthermore, that the radiation field consists of n(< 00) modes, represented 
by creation and destruction operators {a^, ai\l — 0, . .,n — 1} in a Fock-Hilbert space H. ra d as 
defined by the standard specifications that (a) these operators satisfy the CCR, 

[a h a* m ] = 5 lm I] [a h a m ] = 0, (2.7) 

and (b) Tirad contains a (vacuum) vector $, that is annihilated by each of the a's and is cyclic 
w.r.t. the algebra of polynomials in the a*'s. 
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The formal generator of the semigroup T rad of the radiation is 

Lrad = Y^ =0 ( iuj i\- a i a ii •] + 2 «W(-) a z - K>i{a*a h .}), (2.8) 

where {.,.} denotes anticommutator, and the frequencies, u>i, and the damping constants, Ki, 
are positive. We refer to [1] for a rigorous definition of L rad . 

The composite (finite) system is simply the coupled system, T,( N \ comprising the matter 
and the radiation. We assume that its algebra of observables, B^ N \ is the tensor product 
A^t&IZ, where 1Z is the *— algebra of polynomials in the a's, a*'s and the Weyl operators. 
Thus, B^ N \ like 1Z, is an algebra of both bounded and unbounded operators in the Hilbert 
space H (N) := C 4N+2 (g)H ra d- We shall identify elements A, R, of A (N \ K, with A®I rad and 
Imat®R, respectively. 

We assume that the matter-radiation coupling is dipolar and is given by the interaction 
Hamiltonian 

Hint ) =T, relN (^ + 4 N) + h.c.), (2.9) 
where we have introduced the so-called radiation field, 4>^ N \ whose value at the site r is 

0f ) = - t (2N + l)- 1 / 2 ^; =o 1 A z a i exp(2rfr/n). (2.10) 

Here the A's are real-valued, iV— independent coupling constants. 

Among the other results contained in [1], one of the most relevant is that the map 

L^ = L^ t + L rad + l [H^ t \.] 

is really the generator of a iV-depending semigroup, T^ N \ regardless of the unbounded nature 
of both L rad and H-^ t \ This is the starting point for a successive analysis, see [1, 2]. 

Now we introduce the HL model, changing a little bit the notations with respect to the 
original paper, [5] , and introducing n modes for the radiation instead of the only one considered 
by HL. The HL hamiltonian for the 2N + 1 atoms and for the n modes of the radiation can be 
written as follows: 

H = H^ + H( R \ (2.11) 

where " S" refers to the system (radiation+matter) and " R" to the reservoir. The hamiltonian 
of the system is 

n-l n-1 
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a n-1 

+ ^==T E E (*i, + a]e-W» + v^e** 1 '"), (2.12) 



j=o ie/jv 

which differs from the one in [5] for the phases introduced in the last two terms, phases 
which are related to the presence of many modes in this hamiltonian with respect to the original 
one. Notice that the presence of (3 means that we are not restricting our model to the rotating 
wave approximation, (RWA). 

The hamiltonian for the reservoir contains two main contributions, one related to the two 
reservoirs of the matter and one to the reservoir of the radiation. We have: 

if (*) = if + J2 H i A \ (2-13) 

iei N 



where 



rt— 1 „ n— 1 



= £ I dku r ,(k) rj (ky rj (k) + v/S£(rJ(^K + r 3 {g 3 )a]), (2.14) 

3=0 J j=0 

and 

2 

H l A) = E / ^Wm s (^)m! i /(^)^ S ,z(^) + V^(^,K^l)^- + ^-cO + V^(^2,/(^2)^,+ + /?.-C.) (2.15) 
s=l J 

Few comments are necessary in order to clarify the formulas above. 

1) first of all we are using the notation: rj(gj) = J dkrj(k)gj(k) and rj(g~j) = / dk rj(k)g •(/;) . 
Here dk is a shortcut notation for dk 3 . 

2) the functions gj and h± t 2 are introduced by HL to regularize the bosonic fields Tj(k) and 

3) we notice that in this model two indipendent reservoirs, m^i{k) and m,2 t i(k), are intro- 
duced for (each atom of) the matter, while only one, rj(k), is used for (each mode of) the 
radiation. This result will be recovered also in our approach. 

4) it should be pointed out that the hamiltonian above is really only one of the possible 
extentions of the HL original one to the n-modes situation, and, in fact, is quite a reasonable 
extension. In particular we are introducing different dispersion laws and different regularizing 
functions hj for each mode of the radiation, while we use the same u and the same h for the 
atoms localized in different lattice sites. This "non-symmetrical" choice is motivated by the 
AS model itself, where we see easily that the free evolution of the matter observables does not 
dipend on the lattice site, while in the form of the generator (2.8) a difference is introduced 
within the different modes of the radiation. The reason for that is, of course, the mean field 
approximation which is being used to deal with the model. 
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5) the coupling constant y/a is written explicitly for later convenience. 

The role of each term of the hamiltonian above is evident. Instead of using the above 
expression for H, where we explicitly consider the effect of the system and the effect of the 
reservoirs, we divide H as a free and an interaction part, in the following way: 

H = Jfo + vQfj, (2.16) 

where 

n— 1 2 . n— 1 „ 

H = uj r J2 a ) a j + A* E + EE / d ku ms (k)ml ;l (k)m S; i(k) + E / dkuj^ky^k^r^k) 
j=o iei N iei N s=i J j=o J 

(2.17) 

and 

^ = EK r j(#>j + r ife') a j) + E [( m L(^l)^ - + ^ C ") + ("4,1 (^H+ + /i.C.)] + 

j=o «e/jv 



\/a! 



n-1 



= E E K+a,e 2 ^ 7/n + (Ti.-ote- 2 ^'/") + 

J- •>— n r.. 



j=0 /G/jv 
n-1 



E E K+gR 27 ^ + ^.-a.-e^'/n). ( 2 .18) 



^(2^+1) j-oie/^ 

The only non trivial commutation relations, which are different from the ones already given in 
(2.1,2.7), are: 

[r J (k),r l (k'y] = 5 J>l 5(k-k'), [m. %l (k),m\, v {k')] = 6 3>8 >5 l>v 5(k - k') (2.19) 



We end this section by introducing the DHL model. The main difference, which is introduced 
to avoid dealing with unbounded operators, consists in the use of a fermionic reservoir for the 
matter, and for this reason the Pauli matrices of both AS and HL are replaced by fermionic 
operators as described in details, for instance, in [7]. 

The idea for introducing these operators is quite simple: since we are considering only 
two- levels atoms (this is the reason why Pauli matrices appear!) a possible description of 
one such atom could consist in using two pairs of independent fermi operators, (&_,&!) which 
annihilates and creates one electron in the lowest energy level with energy E_, and (b + ,b\_) 
which annihilates and creates one electron in the upper energy level with energy E + . If 
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we restrict the Hilbert space of the single atom to the states in which exactly one electron is 
present, in the lower or in the upper level, it is clear that b\b_ behaves like cr+, that is when 
it acts on a vector with one electron in the lowest state (fel^o being the state with no 
electrons), it returns a state with an electron in the upper level (fe+^o), and so on. Moreover 
b\_b + — has eigenvectors 6±^o with eigenvalues ±1, so that it can be identified with a z . 
Going back to the finite system we put 

a +>l = 6^,6-,,, = &!,,&+,,, a z , t = 6^6+,, - (2.20) 

where I e In- The only non trivial anti-commutation relations for operators localized at the 
same lattice site are: 

{b±,iAi} = 1 - ( 2 - 21 ) 
Moreover, see [7], two such operators commute if they are localized at different lattice site. 
For instance we have [b±j, &± iS ] = if I ^ s. 

Since the number of the atomic operators is now doubled with respect to the HL model, it 
is not surprising that also the number of the matter reservoir operators is doubled as well: from 
2 x (2N + 1) we get 4 x (2N + 1) operators, each one coupled with a b^ 1 operator. On the 
other hand, the part of the radiation is not modified passing from the HL to the DHL model. 
Let us write the hamiltonian for the open system in the form which is more convenient for us 
and using A instead of y/a. We have 

H = H + XHj, (2.22) 

where 

n— 1 n— 1 



It -L I I _L p 

H = u R Y, a)aj + // E ( b +,i b +,i ~ 6 -,i 6 -, + E / dku^ky^r^k) + 

3=0 /G/jv 3=0 J 

+ E E / dke(k)(Bl(k)B.j(k) + Cl{k)CsAk)) (2.23) 



and 

Hi = eVJ(s> + r^a)) + A£ (tf^fc-,, + h.c.) + 

3=0 

+ E E [bU B «(9B.) + CsAgcs)) + (<,(<te.) + 0//rJ)/w. (2.24) 



1 We use here a;" to indicate one of the two possibilities: x or x^ , x being a generic operator of the physical 
system. 
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Here gs± and gc± are real function. 

The commutation rules for the radiation operators (system and reservoir) coincide with the 
ones of the HL model. For what concerns the matter operators (system and reservoirs) the first 
remark is that any two operators localized at different lattice sites commutes, as well as any 
operator of the radiation with any observable of the matter. As for operators localized at the 
same lattice site, the only non trivial anticommutators are 

{B±,i(k), fit ,(*/)} = {C ±>l (k), Ci tl (k')} = 5(k - k'), (2.25) 

together with the (2.21), while all the others are zero. Finally, to clarify the different roles 
between the B and the C fields it is enough to consider their action on the ground state of the 
reservoir <p : 

rj (k)<p = B± tl (k)<p = Clj(k)<p = 0. (2.26) 

These equations, together with what has been discussed, for instance, in [7], show that B is 
responsible for the dissipation, while C is the pump. Again in reference [7] it is discussed which 
kind of approximations, other than using a fermionic reservoir, are introduced to move from 
"real life" to the DHL model. Here we mention only a few: the atom is considered as a two 
level system; only n modes of radiation are considered (n=l in the original model, [5]); the 
electromagnetic interaction is written in the dipolar approximation and within the RWA; the 
model is mean field, etc.. 

It is worth remarking that since all the contributions in H above are quadratic in the 
various creation and annihilation operators, they all commute among them. This fact will be 
used in the computation of the SL of this model. 

Ill Alli-Sewell versus Hepp-Lieb 

We begin this section with a pedagogical note on the single-mode single-atom version of the 
AS model. This will be useful in order to show that two reservoirs must be introduced to 
deal properly with the matter. After that we will consider the full AS model and we will 
show that the hamiltonian which produces the AS generator after considering its SL is nothing 
but the HL hamiltonian in the RWA. We will conclude this section proving that adding the 
counter-rotating term (the one proportional to (3 in (2.12)) does not affect this result, since its 
contribution disappear rigorously after the SL. 

The starting point is given by the set of equations (2.3)-(2.10) restricted to n — 1 and 
iV = 0, which means only one mode of radiation and a single atom. With this choice the phases 
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in (f)\ N ^ disappear so that the interaction hamiltonian (2.9) reduces to 

Hint = i{<y~a ] - h.c), (3.1) 

and the total generator is L = L mat + L rad + i[H int , .}. 

Let us suppose that the atom is coupled not only to the radiation by means of H int , but 
also to a bosonic background m(k) with the easiest possible dipolar interaction: 

H M m = cr + m(h) + h.c. (3.2) 

Of course this background must have a free dynamics and the natural choice is 

H , m = J dkuj m (k)m ] (k)m(k). (3.3) 

For what concerns the radiation background the situation is completely analogous: 

H 0tr = J dku r {k)r\k)r{k), H R , r . = ar^(g) + h.c. (3.4) 

are respectively the free hamiltonian and the radiation-reservoir interaction. We take the 
complete hamiltonian as simply the sum of all these contributions, with the coupling constant 
A introduced as below: 

H = H + XHj = {fia z + u R a^a + J dku m {k)m) (k)m(k) + J dku r {k)r\k)r(k)} + 
+\{(ar j (g) + h.c.) + (a + m(h) + h.c.) + Ai((7_a f - h.c)}. (3.5) 

Taking the SL of this model simply means, first of all, considering the free evolution of the 
interaction hamiltonian, Hj(t) = e tHot Hje~ tHot , see Appendix and reference [6]. It is a simple 
computation to obtain that, if 

lu r = 2/i, (3.6) 

Hr(t) = (ar t (ffe i(w *- Wfl)t ) + h.c.) + {a + m{he ii2 ^ m)t ) + h.c) + Ai(<7_a t - h.c). (3.7) 

In this case the SL produces, see Appendix, the following effective time-depending interaction 
hamiltonian: 

ijj s °(t) = (arj(t) + h.c.) + {cr + m h {t) + h.c.) + i{a_a ] - h.c), (3.8) 

where the dependence on A disappears and the operators r g (t), m h (t) and their hermitian 
conjugates satisfy the following commutation relations for t > t', 

Mf),rJ(f)] = T {9) 5(t - t>), [m h (t),m[(t')} = T W 5(t - f). (3.9) 
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Here we have defined the following complex quantities: 

= f° dr [ dk\g(k)\ 2 e- l ^^-^ T , T W = f° dr f dk\h(k)\ 2 e^ 2 ^™^ . (3.10) 

We want to stress that the restriction t > t' does not prevent to deduce the commutation rules 
(3.13) below, which are the main ingredient to compute the SL. However, the extension to t < t' 
can be easily obtained as discussed in [6]. It is clear now what should be the main analytical 
requirement for the regularizing functions h and g: they must be so that the integrals above 
exist finite! 

In order to obtain the generator of the model we introduce the wave operator XJ t (in the 
interaction representation) which satisfy the following operator differential equation: 

d t U t = -iH\ sl \t)U u with U Q = 1. (3.11) 

In [6], and reference therein, it is proven that for a large class of quantum mechanical models, 
the equation above can be obtained as a suitable limit of differential equations for a A-depending 
wave operator. Analogously, r g (t) and m h (t) can be considered as the limit (in the sense of 
the correlators) of the rescaled operators ±r(ge _i(Wr ~ Wii) * /A2 ) and jm( y he i ^ 2 ^ UJm ^ t/x2 ). It is not 
surprising, therefore, that not only the operators but also the vectors of the Hilbert space of 
the theory are affected by the limiting procedure A — > 0. In particular, the vacuum rj for the 
operators r g and m h , mh(t)rj = r g (t)r) = 0, does not coincide with the vacuum ip for m{k) 
and r(k), r(k)ipo = m(k)ipo = 0, see [6] for more details. 

Equation (3.11) above can be rewritten in the more convenient form 

U t = l-i f Hf(t')U t ,dt', (3.12) 
Jo 

which is used, together with the time consecutive principle, see Appendix and [6], and with 
equation ([?]), to obtain the following useful commutation rules 

[r g (t), U t ] = -iF^aUt, [m h (t), U t ] = -iT [h) a_U t . (3.13) 

If we define the flow of a given observabe X of the system as jt(X) = Xj\XXJ u the generator 
is simply obtained by considering the expectation value of d t jt{X) on a vector state 77^ = 
rjo <8>£, where £ is a generic state of the system. Using formulas (3.11,3.13) and their hermitian 
conjugates, together with the properties of the vacuum i] , the expression for the generator 
follows by identifying L in the equation < d t j t (X) > (o=< j t (L(X)) > (o- The result is 

L(X) =L 1 (X)+L 2 (X) + L 3 (X), 

L X {X) = T { ° ] [a\X]a-T i9) a^a,X}, L 2 (X) = T W [a + , X]a_ -T W a + [a_, X], 

L 3 (X) = i 2 [a_a ] - a + a, X] (3.14) 
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It is evident that both L\ and L 3 can be rewritten in the same form of the radiation and 
interaction terms of the AS generator but this is not so, in general, for L 2 which has the form 
of the AS radiation generator only if the pumping parameter 77 is equal to —1. 

This is not very satisfactory and, how we will show in the following, is a consequence of 
having introduced a single reservoir for the atom. We will show that the existence of a second 
reservoir allows for the removal of the constraint rj = — 1 we have obtained in the simplified 
model above. 

With all of this in mind it is not difficult to produce an hamiltonian which should produce 
the full AS generator for the physical system with 2N + 1 atoms and n modes of radiation. 
With respect to the one discussed above, it is enough to double the number of reservoirs for the 
matter and to sum over / G In for the matter and over j — 0, 1, n — 1 for the modes. The 
resulting hamiltonian is therefore necessarely very close to the HL one: 

H = H + \H h (3.15) 

with 

rt— 1 2 „ n— 1 « 

j=o ' iei N iei N s=i 3=0 J 

(3.16) 

and 

n-1 

H i = J2( r }(9j)aj + r j(9j)a]) + [("^(^iH- + h - c ) + («4#2K+ + h.c.)] + 
3=0 /e/jv 

+AE(#V+ + ^.), (3-17) 
iei N 

where the radiation field has been introduced in (2.10). It is clear that, but for the RWA 
which we are assuming here, there are not many other differences between this hamiltonian 
and the one in (2.11)-(2.15). It is worth mentioning that A appears both as an overall coupling 
constant, see (3.15), and as a multiplying factor of J^ieiNi^i^ 171 ^ + h.c.) and plays the same 
role as y/a in the HL hamiltonian. As for the commutation rules they are quite natural: but 
for the spin operators, which satisfy their own algebra, all the others operators satisfy the CCR 
and commute whenever they refer to different subsystems. In particular, for instance, all the 
m[ tl (k) commute with all the m\^(k'), for all k, k' and /, I'. 

The procedure to obtain the generator is the same as before: we first compute Hj(t) = 
e iH tjj ie -iH t^ enters in the differential equation for the wave operator. Taking the limit 
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A — > of the mean value in the vector state defined by ipffi = tp <g> £ of the first non trivial 
approximation of the rescaled version of U t we deduce the form of an effective hamiltonian, 
Hj Sl \t), which is simply 

n— 1 

Hi sl \t) = E(%4W+^- c -)+EK+^iX^)+^-c-)+EK-^(0+^-c.)+E (4 N) <n, + +h.c). 

j=o iei N leiN ie-fjv 

(3.18) 

Again, we are assuming that ujr = 2/x, which is crucial in order not to have a time dependence 
in the last term of Hf\t) in (3.18). 

The only non trivial commutation rules for t > t' for the new operators are: 

[r,j(*),ry«01=r?i^(*-«0, 

Ki(*),™U*')] = rL h %,<5(* - 0, (3-i9) 

where we have defined the following complex quantities: 

= f dr f dk\ gj (k)\ 2 e l ^ k) -^ )T , = f dr f dk^k^e^^^ , 

J — oo J J —oo J 

T (M) = f° dr f dk\h 2 (k)\ 2 e l ^ m ^ +w ^ T . (3.20) 

J -co J 

The above commutators, given for t > t', are sufficient to compute the commutation relations 
between the fields of the reservoir and the wave operator U t — 1 — % J Hj Sl \t')U t >dt' , as it is 
obtained after the SL. We get 

[r g j(t), U t ] = -iT {g ] 3 a 3 U u [m hl (t), U t \ = -iT^a^Ut, [m 2i/ (t), U t \ = -iT™ <r li+ U t . (3.21) 

The expression for the generator can be obtained as for the N = 0, n = 1 model described 
before, that is computing the mean value < dtjAX) > Here, as before, jt(X) is the flux 
of the system obserable X, jt(X) = U}XU t , and i] is the vacuum of the operators r g j(t) and 
m s j(t), s = 1,2. The computation gives the following result, which slightly generalize the one 
in (3.14): 

L(X) = L 1 (X) + L 2 (X) + L 3 (X), 

ra-1 

MX) = E( r -iH*K- -r^-a]K-,x]), 

j=0 

L2(X) = Y: (T^\a + , l ,X]a^-T {hl) a + ,[a^ l ,X]+T^\a_, l ^ 

L 3 (X) = i Y, + h.c.),X}. (3.22) 
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It is not difficult to compare this generator with the one proposed by AS, see formulas ((2. 3), (2. 10)), 
and the conclusion is that the two generators are exactly the same provided that the following 
equalities are satisfied: 

31^ = ^, 1tr"j /,> sR(T (hl) + T {h2) ) = 7l , %(T (hl) -T {h2) )=e 

KTL fcl) = ^72(1-17), ^ = ^ 2 (1+ V ). (3.23) 

Here the lhs all contain variables of the hamiltonian model while the rhs are related to the AS 
generator. Due to the fact that the hamiltonian in (3.15)-(3.17) essentially coincides with the 
one in (2.11)-(2.15) with f3 — 0, that is in the RWA, we can conclude that if we start with the 
HL hamiltonian, choosing the regularizing functions in such a way that the equalities (3.23) 
are satisfied, the SL produces a generator of the model which is exactly the one proposed in 
[1, 2], with the only minor constraint 72 = 271, which is a direct consequence of (3.23). 

From a physical point of view the implications of this result are quite interesting: the original 
model, [5], was not (easily) solvable and for this reason a certain number of approximations 
were introduced. Among these, the crucial ones are the replacement of the original reservoir 
with what HL call a singular reservoir which, moreover, is made of fermions. Under these 
assumptions the model can be discussed in some details, and this was done in [8], where the 
thermodynamic limit for the intensive and the fluctuation observables was discussed. What 
we have shown here is that all these approximations can be avoided using another kind of 
perturbative approach, that is the one provided by the SL. The resulting model is exactly the 
one proposed and studied in [1, 2]. The role of the singular reservoir, or the need for a fermionic 
reservoir, is therefore not crucial and can be avoided. However, we will consider the DHL model 
in the next section in order to complete our analysis. 

We conclude this section with a remark concerning the role of the RWA and its relation with 
the SL. In particular, this is a very good approximation after the SL is taken. To show why, we 
first notice that adding a counter-rotating term (extending the one in (2.18)) to the interaction 
hamiltonian Hj in (3.17), considering the same coupling constant for both the rotating and 
the counter-rotating term ((3 = a), simply means to add to Hj in (3.17) a contribution like 
A J2i£i N (4 , i N ^ <J i,- + h.c). While the rotating term, if 2/j, = ujr, does not evolve freely, the free 
time evolution of this other term is not trivial. However, the differences with respect to the 
previous situation all disappear rigorously after the SL, because these extra contributions to 
the mean value of the wave operator go to zero with A, so that at the end the expression for the 
generator is unchanged. This allows us to conclude that the full HL hamiltonian is equivalent 
to the AS generator, where the equivalence relation is provided by the SL. 
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We want to end this section with a final remark concerning the different number of phase 
transitions in the two different situations, 1 for the HL and 2 for the AS model. In view of 
the above conclusion, we can guess that the SL procedure produces some loss of information 
and, as a consequence, a difference between the original and the approximated system. This 
is not so surprising since, thought being a powerful tool, nevertheless the SL is nothing but a 
perturbative method! 

IV The SL of the DHL model 

In this section we consider the SL of the DHL model as introduced in Section 2. In particular we 
find the expression of the generator and we show that, under some conditions on the quantities 
defining the model, the equations of motion do not differ from the ones in AS. The free evolved 
interaction hamiltonian Hj in (2.24) is, 

n— 1 

H T (t) = e iHot H ie - iHot = £(%rj(^--^)*) + h . c .) + A £ (4*%^ + h.c.) + 

j=o N 

+ E [bUB^(g B+ e^) + C + ^g c+ e^)) + (B^g^e-*^) + C^^+e-^ + 
i&i N 

+bl l (B^ l (g B _e-^) + C^{gc-e- lt ^)) + (<^-^ +e) ) + C^gc-e^)^}. (4.1) 

Following the usual strategy discussed in the Appendix and in [6], we conclude that (the rescaled 
version of) the wave operator U\(t) = 1 — i\ / * Hj(t')U\(t')dt' converges for A — > to another 
operator, which we still call the wave operator, satisfying the equation 

U t = l-i / Hf s \t')U t ,dt', or, equivalently d t U t = -iHf s \t)U t , with U = 1. (4.2) 
Jo 

Here Hf s) (t) is an effective time dependent hamiltonian defined as 

H? s) (t) =£WL(*) + h - c -) + Y, lelN (<Pl N KiKi + h.c.) + £ [b^(Mt) + 7+ ,^)) + 
o=o iei N 

+(PUt) + ilMKi + bUP-At) + -r-M) + (P-M + iUW-A- ( 4 - 3 ) 

The operators of the reservoir which appear in 

H (ls) 

are the stochastic limit of the original 
(rescaled) time evoluted operators of the reservoir and satisfy (anti-)commutation relations 
which are related to those of the original ones. In particular, after the SL, any two operators of 
the matter (system and reservoirs) localized at different lattice sites commutes, as well as any 
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operator of the radiation with any observable of the matter. As for operators localized at the 
same lattice site, the only non trivial anticommutators are 

{/MO, 4/0} = 6(t - t')ri B± \ { 1± MlUt')} = 5(t - t')T {c± \ (4.4) 

which should be added to 

[r a At),r B At)] = 8jsS(t-l0r% (4.5) 

In all these formulas the time ordering t > t' has to be understood and the following quantities 
are defined: 

r^- = [° dr [ dk\ gj {k)\ 2 e i{ - w rAk)-"R)r ? r (f±) = f° dT f dk(g B± (k)) 2 e l « k) ^ )T , 
T (C±) = J° dr J dk(g c ±(k)) 2 e- l{t{k)T>l)T . (4.6) 
We call now rj the vacuum of these limiting operators. We have 

r g ,Mvo = P±,i(t)m = ^f±,i(t)Vo = 0. (4.7) 

Paying a little attention to the fact that here commutators and anti-commutators simultane- 
ously appear, we can compute the commutators between the operators r 9! j(t)^, 7± j(£)> P± i{t) 
and the wave operator U t by making use of ((4. 4), (4. 5)). We give here only those commutation 
rules which are used in the computation of the generator: 

KM U t ] = -iT^ajUt, \P±M U t ] = -iT (B±) b ±il U t , frlM U t ] = tT^b^Ut. (4.8) 

We are now ready to compute the generator following exactly the same strategy as in the 
previous section and in Appendix. However, we must remark that, due to the presence of 
fermionic operators, to make the computation simpler, we will focus on system operators X 
which are quadratic (or quartic,...) in the matter operators localized in a given lattice site (e.g. 
X = b+jb'Lj, X = b+jfi+j^jb'Lj,...) in such a way to ensure the commutativity between X 
and any of the matter operators entering in the operator Hj ls \t), (4.3). We get for the full 
generator the following sum of three different contributions: 

L(X) = L 1 (X) + L 2 (X) + L 3 (X), 

n—l 

Li(X) = £(r^hU]a, -T%a}\ aj ,X]), 

3=0 
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L2(x) = £ (r^ib^xjb+j - r { _ B+) bl,[b + , h x] + rL° - r^b^b^x] + 

+Y^ B -\bl p X]b^-Y {B -hl l [b^X] + Y^ 

L 3 (X)^i^U < i N Kib-,i + h.c.),X}. (4.9) 
lei N 

We see that the first and the last terms exactly coincide with the analogous contributions of the 
AS generator, but for a purely formal difference which is due to the different matter variables 
which are used in the two models. The second contribution, on the other hand, cannot be 
easily compared with the free AS matter generator. What is convenient, and sufficient, to get 
full insight about L 2 , is to compute its action on a basis of the local algebra, that is on b\ fi^ i 
(= a +! i) and on b^ + l b +t i — 6_ (= &z,i), an the others being trivial or an easy consequence of 
these ones. It is not hard to find the result: 

Hbi^,) = + r<f-> + rL c+) + r<?->] - 

-^[r^-r^-r^ + r^]), 

Little - bljb-j) = 2(-bl il b + , l (^- +] + + ^- +) + 

+bl^ l mv (B - ) + rL c_) ) - 3?rL c_) ). (4.10) 

The equation for a + j is recovered without any problem, modulo some identification (K[ll B+) + 
r ( i B_ ' ) + r < L 7+ - ) + r^ - * 1 ] = 7 1; ...), while to recover the equation for a z j it is necessary to choose 
properly the regularizing functions which define the different T_. In particular we need to have 
the following equality fulfilled: 

£( r (f+) + r {c+) ) = K(r (B_) + T {c ~ ] ). (4.11) 

Under this condition we can conclude that the SL of the DHL model produces the same dif- 
ferential equations as the AS generator, as for the HL model. It is also easy to check that, as 
a consequence of our approach, we must have 71 = 72 in the generator we obtain. Of course 
this result is not surprising since already in the HL paper, [5], the fact that the two models are 
quite close (under some aspects) was pointed out. Here we have learned also that the SL of 
both these models, at least under some conditions, give rise to the same dynamical behaviour. 



V Outcome and Future Projects 

We have proved that the relation between the HL and the AS model, whose existence is claimed 
in [1], is provided by the SL. This result is quite interesting since it shows that the approx- 
imations introduced by HL in their paper [5], in particular the use of the fermionic reservoir 
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for the matter which produces the DHL model, together with the so called singular reservoir 
approximation, can be avoided by using the original HL model with no approximation, taking 
its SL and finally using the results in [1, 2] to analyze, e.g., the thermodynamical limit of the 
model. It is interesting to remark that while in the AS model two phase transitions occur, in 
the HL model we only have one. This could be a consequence of the SL procedure, which is 
nothing but a perturbative approach simplifying the study of the quantum dynamics, so that 
some of the original features of the model can be lost after the approximation. 

We want to conclude this paper by remarking that this is not the first time the HL model is 
associated to a dissipative system, as in the AS formulation. A similar strategy was discussed 
by Gorini and Kossakowski already in 1976, [9]. It would be interesting to study their generator 
again in the connection with the SL to see if any relation between their generator and the HL 
original hamiltonian appears. 
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A Appendix: Few results on the stochastic limit 

In this Appendix we will briefly summarize some of the basic facts and properties concerning 
the SL which are used all throughout the paper. We refer to [6] and reference therein for more 
details. 

Given an open system S+R we write its hamiltonian H as the sum of two contributions, the 
free part H and the interaction XHi. Here A is a (sort of) coupling constant, H contains the 
free evolution of both the system and the reservoir, while Hi contains the interaction between 
the system and the reservoir and, for composite systems, amoung the different buildind blocks of 
the whole physical system. Working in the interaction picture, we define Hi(t) = e' lH ° t H I e~' lHot 
and the so called wave operator U\(t) which satisfies the following differential equation 

d t U x (t) = -i\H I (t)U x (t), (A.l) 

with the initial condition U\(0) = 1. Using the van-Hove rescaling t — > see [7, 6] for instance, 
we can rewrite the same equation in a form which is more convenient for our perturbative 
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approach, that is 

W X (±) = -i^(l)C/ A (A), (A.2) 
with the same initial condition as before. Its integral counterpart is 

U x {±) = l-if o H I {l)U x {l)dl?, (A.3) 

which is the starting point for a perturbative expansion, which works in the following way: let 
tpo be the ground state of the reservoir and £ a generic vector of the system. Then we put 
ip^ — tpo <S> £. We want to compute the limit, for A going to 0, of the first non trivial order of 
the mean value of the perturbative expansion of U\(t/\ 2 ) above in ip^ \ that is the limit of 

h(t) = (-i) 2 jT* dhj^dh <H I {^)H I {^) >^), (A.4) 

for A — > 0. Under some regularity conditions on the functions which are used to smear out 
the (typically) bosonic fields of the reservoir, this limit is shown to exist for many relevant 
physical models, see [6] and [10] for a recent application to many body theory. At this stage 
all the complex quantities like the various T_ we have introduced in the main body of this 
paper appear. We call I(t) the limit \im\_> I\(t). In the same sense of the convergence of the 
(rescaled) wave operator U\{j^) (the convergence in the sense of correlators), it is possible to 
check that also the (rescaled) reservoir operators converge and define new operators which do 
not satisfy canonical commutation relations but a modified version of these. Moreover, these 
limiting operators depend explicitly on time and they live in a Hilbert space which is different 
from the original one. In particular, they annihilate a vacuum vector, i] , which is no longer 
the original one, <po- 

It is not difficult to write down, as we have done several times in this paper, the form of a 
time dependent self-adjoint operator Hf s \t), which depends on the system operators and on 
the limiting operators of the reservoir, such that the first non trivial order of the mean value 
of the expansion of U t — 1 — i Jq Hj ls \t')U t >dt' on the state r]^ = 7] <g> £ coincides with I(t). 
The operator U t defined by this integral equation is called again the wave operator. 

The form of the generator follows now from an operation of normal ordering. More in 
details, we start defining the flux of an observable of the system X = X <g> l res , where l res is the 
identity of the reservoir, as j t (X) = U ] t XU t . Then, using the equation of motion for U t and Lf}, 
we find that d t jt(X) = iU}[Hf s \t) , X]U t . In order to compute the mean value of this equation 
on the state r]^\ so to get rid of the reservoir operators, it is convenient to compute first the 
commutation relations between U t and the limiting operators of the reservoir. At this stage the 
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so called time consecutive principle is used in a very heavy way to simplify the computation. 
This principle, which has been checked for many classes of physical models, and certainly holds 
in our case where all the interactions are dipolar, states that, if (3(t) is any of these limiting 
operators of the reservoir, then 

\P(t),U v \ = 0, for all t > t'. (A.5) 

Using this principle and recalling that i] is annihiled by the limiting annihilation operators 
of the reservoir, it is now a technical exercise to compute < d t jt(X) > ^ and, by means of 
the equation < d t j t (X) > to=< j t (L(X)) > («, to identify the form of the generator of the 
physical system. 
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